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Abstract 

Both the mass (just below the D*^D^ threshold) and the likely quantum numbers {J^'-'' = 
of the X(3872) suggest that it is either a weakly-bound hadronic "molecule" (X(3872 ~ 
l/\/2[D*^D^ + D*^D^]) or a virtual state of charmed mesons. Assuming the X(3872) is a weakly- 
bound molecule, the scattering of neutral D and D* mesons off the X(3872) can be predicted from 
the X(3872) binding energy. We calculate the phase shifts and cross section for scattering of D'^ 
and D*^ mesons and their antiparticles off the X(3872) in an effective field theory for short-range 
interactions. This provides another example of a three-body process, along with those in nuclear 
and atomic systems, that displays universal properties. It may be possible to extract the scattering 
within the final state interactions of Be decays and/or other LHC events. 
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I. INTRODUCTION 



In recent years many new and possibly exotic charmonium states have been observed 
at the B-factories at SLAC [2], at KEK ^ in Japan, and at the CESR coUider at Cornell 
[3|. This has revived the field of charmonium spectroscopy 0, 0, 0, 0, B [9|, which will be 
an important part of PANDA at the FAIR facility [lo|. Because several of the new states 
exist very close to scattering thresholds, it is useful to interpret them as hadronic molecules, 
a concept introduced in Refs. (ill . IT^ . 13 1 well before these most recent experiments. A 



particularly interesting example is the X(3872), discovered by the Belle collaboration [14 



in —* K^-K^Ti J/ip decays and quickly confirmed by CDF 15|], DO [16|, and BaBar 
[17(1 . The state has likely quantum numbers J^^ = and is very close to the 0*^0^^ 
threshold.^ As a consequence, the X(3872) has a resonant S-wave coupling to the D*^D^ 
system. Early examples of discussions of the possible molecular nature of the X(3872) can 



be found in Refs. [18j, [19|, l20| . An extensive program (see a status report in Ref. [2l|]) provides 
predictions for its decay modes based on the assumption that it is a D*^D^ molecule with 
even C-parity: 
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(1) 



The measured mass and width of the X(3872) differ significantly in the J/ipTT^n' and 
£)*ojjo jgQay channels. This effect can be understood from a line shape analysis which shows 
that the true mass and width of the X(3872) are measured in the Jlijj7r^7r~ channel because 



the D*^D^ channel is contaminated by a threshold enhancement ^ ^ 

latest measurements in the J /ipTi^Ti~ channel [IgI, 27, 23], the mass of the X(3872) is 2^: 




24,l25|. Usingthe 



mx = (3871.55 ±0.20) MeV, 
which corresponds to an energy relative to the D*^D^ threshold ^1 of 



(2) 



(3) 

0.26 

and 

3ll . l32j). The X(3872) is also very narrow, with a width smaller than 



Ex = (-0.26 ±0.41) MeV. 

The central value corresponds to a {D*'^D^)^ bound state with binding energy Bx 
MeV (but a virtual state cannot be excluded from the current data in the J/ipn^ 
D*^D^ channels |P 
2.3 MeV. 



TT 



Because the X(3872) is so close to the 



threshold, it has universal low-energy 
Close to threshold, the coupling to 
threshold is about 8 MeV higher in 



properties that depend only on its binding energy [33 . 
charged D mesons can be neglected because the D*^D~ 
energy. Therefore, the properties of the X(3872) can be described in a universal EFT with 
contact interactions only. This EFT is a pionless EFT because the pion degrees of freedom 
are not dynamical near threshold; they are integrated out and all effective interactions are 



short range. This EFT is widely used in low-energy nuclear physics |34l. l35l. l36l. l37l| . For an 



EFT of the X(3872) including explicit pions, see Refs. I 

as a (D*°Z)°)+ molecule in the pionless EFT was initiated by Braaten and Kusunoki [40i]. A 



Q. The study of the X(3872) 



number of predictions for production amplitudes [42], decays [43|, formation [41j, and line 



^ Note, however, that J^'^ = 2 + cannot be experimentally excluded at present. 
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shapes [22|, 23 1 within this framework followed. The interactions of the X(3872) with other 
hadrons, however, are not known. 

In this paper, we extend these studies to three-body processes in the pionless EFT. Based 
on the assumption that the X(3872) is an S-wave (Z}*''D°)+ molecule, we provide model- 
independent predictions for the scattering of and D*^ mesons and their antiparticles 
off the X(3872) in the pionless EFT. We will refer to these reactions collectively as D^*^X 
scattering. This reaction may contribute to the final state interaction in decays of mesons 
into D and D* mesons, in rare events in BB production where one of the -B's decays into 
an X and the other one into a or D* meson, and in prompt events at colliders. In the 
next section we will provide the EFT for D^*^X scattering. In section IIIII we will present 
our results and discuss possible scenarios for observing this process. 



II. FORMALISM AND CALCULATION 

In this section we set up the EFT for D^*^X scattering, derive the integral equation 
for the scattering amplitude, and provide an expression for the total cross section. The 
formalism can be taken over from the pionless theory in low-energy nuclear physics, but we 
briefly describe the issues relevant for D^*^X scattering. For a more detailed discussion and a 



bibliography of the original work, see the reviews of Refs. [3J, |35|, 136|, [371]. For the derivation 
of the three-body equations, it is convenient to introduce a non-dynamical auxiliary field X 
for the X(3872). The EFT is organized in an expansion around the non-trivial fixed point 
of the coupling between and Z)°* mesons corresponding to infinite scattering length or, 
equivalently, the X(3872) being a threshold bound state. The binding momentum, 7, of the 
X(3872) is y/2fixBx, with the reduced mass nx = mj:,om£,*o / {m^fO + m£,*o). So an infinite 
scattering length corresponds to the 7 = limit. The EFT expansion is then in powers of 
7/A6 and k/ Kb where k is the typical momentum exchange and A;, is the breakdown scale 
of the pionless EFT. We will estimate A;, from one-pion exchange in the discussion of the 
errors below. To leading order in this expansion, the effective Lagrangian for the interaction 
of the X(3872) with neutral D and D* mesons can be written as: 



[X^^do^Pd^o + V-d^oVdo) + H.c. ) + ..., (4) 



j=DO,D*o^DO,D* 

9_ 

V2 



where H.c. denotes the Hermitian conjugate and the dots indicate higher order terms with 
more derivatives and/or fields. The terms with more derivatives are suppressed at low 



energies. As shown in [40|, there is no Efimov effect [4J] in this system, so three-body 
terms will not contribute up to next-to-next-to-leading order (N2L0) in the expansion in 
7/A;,. Four- and higher-body forces do not contribute in three-body processes. The only 
corrections up to N2L0 are effective range contributions and their inclusion is in principle 
straightforward 45, 46, 13]. However, the effective range for the X(3872) is not known. We 



will therefore restrict our calculation to leading order only, and estimate the size of higher 
order corrections. The D^, D*^, Z)°, and D*^ mesons are treated as distinguishable particles; 
charge conjugation invariance yields m/jo = m^o and 1710*0 = rrifj^o. 

The parameters A and g in Eq. (jl]) are not independent; only the combination g'^/A 
enters into physical observables. Since the theory is nonrelativistic, all particles propagate 
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FIG. 1: Dressing of the bare X propagator (double dashed hne) by and D*^ meson loops (solid 
and dashed lines). 



forward in time and the tadpoles vanish. The propagator for the D^*^ mesons is 

tS,{po, p) = , ' , J = D\ D', D*' , 

Po ~ p /[2mj) + ze 



(5) 



where = p^. The X propagator is more comphcated because of its couphng to two-meson 
states. The bare X propagator is constant, iDx,o{po,p) = i/^, but the full propagator 
is dressed by and D*^ meson loops to all orders (see Fig. [1]). The bare and full X 
propagators are indicated by the double dashed and double lines, respectively. The D 
mesons are indicated by the solid {D^ and D^) and dashed (Z)*° and D*^) lines. Note that 
each loop receives contributions of two combinations of D^*^ mesons: D^D*'^ and D^D*^ . 
Summing the resulting geometric series leads to the full X propagator: 



iDxipo,p) = iDxfi{po,p)[l - Dxfi{po,p)^{po,p)] ^ 



(6) 



where S(p) is the self energy of the X. Using the reduced mass of the and D*^ mesons 
/ix and their total mass Mx = frijjo + m£)*o, the self energy can be written 



S(po, P) 



2fix9' 

4:71 



(2^ 



M 



-, -1 



p ■ q 



le 



X 



It 



-A + 0(1/A) 



TT 



(7) 



where the ultraviolet divergence was regulated with a momentum cutoff A. Substituting 
this expression into Eq. ([6]) and dropping terms that vanish as A — >■ cxo, we obtain the full 
X propagator: 



iDx{po,v) 



-7 + 



2ixxg^ 

where we have matched the bound state pole position 

1 47r A 2 



M 



It 



X 



(8) 



7 



A. 



TC 



(9) 



to the binding momentum of the X(3872): 7 = ^/2fIxBx and a is the D^D*^ scattering 
length. The combination of bare coupling constants (^^/A must depend on the cutoff A as 
prescribed by Eq. ([9]) since a and 7 are physical quantities. 

Using the full X propagator we may calculate the scattering of a or D*^ meson (or 
their antiparticles) off the X(3872). Because of their different masses, the scattering of a 

versus a, D*^ will lead to a different scattering amplitude and cross section even though 
the interaction strength (^^/A is the same. The scattering amplitude is the solution of the 
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FIG. 2: Integral equation for scattering of a particle S (single line) off the X(3872) (double line). 
The dashed line indicates the particle S complementary to S as explained in the text. 



integral equation shown in Fig. [21 The X and the scattered meson, denoted by S, are 
represented by a double line and a single line, respectively. If an S particle is scattered, a 
second complementary particle type, S (represented by a dashed line in Fig. [2]) appears in 
the loops. The masses of the 5* and S particles are different. For a given scattered particle 
S, the corresponding particle S complementary to S can be read off the flavor wave function 
of the X(3872) in Eq. ([1]). For example, ii S = then S = D*^. 

We now formulate the scattering problem in the center-of-mass frame of the SX system. 
With k the relative momentum of 5* and X, the total energy is 

E=-^-Bx, (10) 

where fisx = ^nsMx / {Mx+ms) is the reduced mass of the SX system and S = , _D°, 
or Z)*°. The resulting integral equation for the off-shell SX scattering amplitude is 

T(k, p) = ^ + fc^ + ■ k - 2/.^^' 



(277)2 J ^ ^-2iix iE~qy{2fisx)) 



It 



X (/ + g2 + ^p-q-2^xi5;) , (11) 

V "^s / 

where k and p are the relative momenta in the incoming and outgoing SX system, respec- 
tively. Performing a partial wave decomposition of T(k, p), 

T(k,p) = ^(2/ + l)T,(A;,p)P,(cosM> (12) 

where Q^p is the angle between k and p and Pi{cos9kp) is a Legendre polynomial, and 
projecting onto the l-th partial wave, we obtain 

m,p) = '-^i^A-lYQi (^Ap' + k'-2f^xE) 
IJ,x 2fixpk \2fixpk 



^1 



^ Jo -7 + ^/^2^^^^[W^^¥JW^x))^^ 

x^{-m(^{p' + q'-2^^xE)] , (13) 
2^xpq K^f^xpq J 
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where 



(14) 



is a Legendre function of the second kind. The integral equation for the S-wave amphtude 
To^kjp) reduces to 



27r7 ms , ( f + k'' + '-j^^pk - 2f,xE 



To{k,p) = ^-^In 



Hx 4:jj,xpk \ p2 + - i^pk - 2fixE 

1 r,, o'm,,) 

vr ■ 



-7 + {E - qy{2fisx)) - ie 



X ^ ' ' " ] (lb) 

^^ixpq \p^ + g2 - it^pq - ■ ^ ^ 

Solutions of the integral equations f|T3|) and f|T5|) can be obtained numerically using standard 
techniques. 

The amplitudes Ti are related to the scattering phase shifts through the relation: 

Ti{k,k) = ^ ^ -. (16) 

jjsx k cot d[k)i — ik 



p' + q' + ^pq-2fixE' 



Using the expression for the differential cross section in terms of the phase shifts: 

(17) 



da 

dn 



> ; rr-n cos 

^kcot5i-ik ^ 
I ' 



we obtain the total cross section for SX scattering: 

TT 



axs{E) = J2^^^±^m{k,k)\\ (18) 
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III. RESULTS AND DISCUSSION 



In this section we present our results for the SX scattering amplitude and the total cross 
section at leading order (S = D°,D*°,D^, or D*^). With the masses of the and D*^ 
mesons fixed, these quantities depend at this order only on the binding momentum 7 of the 
X(3872) (or, equivalently, the D*^D^ scattering length a = I/7). Our results are given in 
units of the scattering length and may be scaled to physical units once a is determined. At 
present the error in the experimental value for Ex in Eq. implies a large error in the 
scattering length. In particular, we obtain the ranges 7 = (0 — 36) MeV and a = (5.5 — co) 
fm with central values 7 = 22 MeV and a = 8.8 fm. 

In Fig. [21 we show our results for the S-wave scattering amplitude fo{k) = l/{k cot 5o(^) ~ 
ik) for the scattering of and D*^ mesons off the A(3872) for center-of-mass momenta 
k from threshold up to O.57. These momenta are still well below the breakup threshold 
oi ks = 1.147 for D^X scattering and ks = 1.1 77 for D^*X scattering. The scattering 
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FIG. 3: S-wave scattering amplitude fo{k) = 1/ {k cot 6o{k) — ik) for scattering of and D*^ 
mesons off tlie X(3872) in units of the scattering length a. The scattering amplitude is identical 
for particles and antiparticles. 



amplitude of a particle is the same as that of its corresponding antiparticle. There is clearly 
an enhancement of the real part of the scattering amplitude at threshold by a factor of 10 
to 17 depending on whether D^X or D*^X scattering is considered. This could lead to an 
enhancement of the scattering cross section by two orders of magnitude compared to the 
already large cross section for D^D*^ scattering. The resulting scattering lengths for D^X 
and D*^X scattering are 

c^D0x = ~9.7a, and a£)*ox = — 16.6 a . (19) 

The corresponding cross sections as a function of the center-of-mass momentum k are 
shown in Fig. HI The difference between the contribution of S-waves (/ = 0) and the full 
cross section (including all partial waves up to / = 6) is negligible for momenta below 7. 
Using the central value of the scattering length a estimated above, we obtain for the scale 
factor a? = 0.78 barn. This factor can become infinite if the X(3872) is directly at threshold, 
while the lower bound from the error in Ex would give a value of 0.3 barn. Even in this 
case the total cross section at threshold will be of the order 300 barns for D'^X scattering 
and 1000 barns for D^*X scattering. 

Next we will estimate the leading corrections to our results arising from the effective 
range of the D^D*^ system. The effective corrections can in principle be calculated up to 



next-to-next-to- leading order in our expansion in a straightforward way 45|, |46|, |47j, but 
since the effective range is not known for this system, an estimate will suffice. 

Naively one expects the breakdown scale of the pionless theory (and the size of the 
effective range) to be set by the pion mass just as in nucleon-nucleon scattering since the 
longest range interaction not explicitly included is the one-pion exchange. The situation 
in the D^D*^ system is potentially more interesting |i38i] because the mass splitting of the 
and D*^ mesons, A = 142 MeV, is almost of the same size as the neutral pion mass. 
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D-X (L=0) 

D*-X(L=0) 

--■ D-X (L<7) 

-- D*-X(L<7) 




0.5 1 
k[1/a] 

FIG. 4: Total cross section for scattering of and D*^ mesons off tlie X(3872) for S-waves (/ = 0) 
and including higher partial waves with Z < 7, in units of the scattering length a. The cross section 
is the same for the scattering of particles as it is for the scattering of antiparticles. 



= 135 MeV. The range of the one-pion exchange interaction is set by the smaller scale /i = 
v/A2 - ml ^ 45 MeV. The mass difference A appears in the propagator of the exchanged 
pion because it carries energy g° = A, leading to the one-pion exchange amplitude for the 
£)0£)*o interaction 



38 



g e ■ qe ■ q 



q2 — /i^ 



(20) 



where g is the D-meson axial transition coupling, f^^ the pion decay constant, e and e* the 
polarization vectors of the incoming and outgoing D^* mesons, and q the three-momentum of 
the exchanged pion. However, the work of Ref. ^] shows that, in part because of the small 
size of the axial coupling, this contribution to NLO effects is very small; contact corrections 
will dominate, bringing us back to an estimate of the effective range of tq ~ I/ttl.^ ~ 1.5 
fm (Af, ~ m7r). Since the leading corrections to our results are of order kvQ and ro7, we 
expect the errors to remain less than 20% even at momenta close to breakup (using the 
central value of the binding energy of the X(3872) in Eq. ([3])). For larger momenta, the 
error is dominated by the kvQ correction and will increase to 35% at momenta of order 45 
MeV. Compared to the errors from effective range corrections the effects from the charged 
D meson channel can safely be neglected; they only enter at much higher momenta since the 
energy difference between the neutral and charged thresholds of about 8 MeV corresponds 
to typical momenta of order 130 MeV or 67. 

To observe the three body interactions described here requires identifying an experimental 
process where, for example, two mesons and one D*^ are produced very near each other 
in space and time. One possibility is provided by the decay of the particle. The B^. 
was discovered through its decays into J/ ip in Run I at CDF (isl . |49| . Particle Data Book 
(2007) averages are: rriB, = (6.286 ± 0.005) GeV and tb, = (0.46 ± 0.07) x lO^^^ s_ Several 
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c ^ c 

FIG. 5: An example of a quark-level B^. decay yielding three charmed/anticharmed quarks in the 
final state 



analyses have been undertaken (see references in [SOj, Chapter 4) to determine the most 
likely mode by which the would decay; the b quark decaying first, the c quark decaying 
first, the two valence quarks annihilating, etc. For access to the three body neutral D^*'^ 
meson interactions, we require that the decay in a mode such as that in Fig. 0, yielding 
three c(c) quarks in the final state. The mass total for the three body D meson system will 
be 5.75 to 5.88 GeV (depending upon whether the third Z) is a or a D*^). Along with the 
additional meson in a P-wave required to balance the Be charge and spin, there is not much 
phase space available. The q in the diagram could be the Cabbibo-favored strange quark 
or the Cabbibo-suppressed {Vcd/Vcs ~ 1/4) down quark. Relative suppression of both decay 
modes is caused by Pauli interference between the spectator c and the c from the (second 
vertex of the) weak decay of the_6 quark. From Ref. [sO]: the estimate for the quark level 
Be cccs decay is about 1.4% |51|]; the detection efficiency for a single is expected to be 
11 to 31%; the production cross section at the LHC (not including feeddowns, which may 
provide an increase of more than a factor of five, but might also be harder to identify) of the 
Be is expected at the 30-60 nb level. At LHCb, the yield will be perhaps 10'' Be events per 
week of running. So the prospect of seeing the three body neutral D^*^ meson interactions 
through Be decay may well be difficult, but is worth investigating. 

Another possibility for observing XD^*^ scattering would be in a BB production event 
where one B decays to an X(3872) and its partner B decays to a D^*^ [52]. Conditions at 
asymmetric 5-factories do not favor the interaction of the decay products [53], but conditions 
at the LHC may. Heavy flavor production at the LHC has received extensive attention 
because of the need to correct for Standard Model background processes in the search for the 



Higgs or new physics. The theoretical and experimental prospects are reviewed in Ref. [54 
Production of the X(3872) under CDF conditions is already dominated by prompt events 
[55] , and heavy flavor production at the LHC is expected to be dominated by gg fusion. The 
cross section for bb at the LHC is about 0.5 mb. Predictions of bb correlations indicate that 
there are events where the opening angle between them may be small (see Fig. 5 on pg. 266 
in Ref. 13]); the gg QQ cross section is dominated by rapidity differences less than one 



[56|. This is important if we expect their decay products to interact. The cross section for 
cc production at the LHC is larger (10 mb), but a minimum of cccc would be required to 
produce an XD^*^ scattering event. The ALICE detector will be sensitive to quarkonia-type 
particles such as the X(3872), as well as D particles, which might be produced promptly, 
while the LHCb detector is optimized to look at B decay products at larger rapidity, where 
an XD^*^ interaction might be more likely. Another process that might yield smaller opening 
angles and an enhanced opportunity for final state interactions would be the production of 
X and D^*^ from bbbb, which is expected at the LHC at a cross section of ~ 500 fb (57| . 
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Nucleus- nucleus collisions in the LHC will also produce X(3872)'s along with associated 
D^*^ mesons. 

The effect of final state interactions involving the X and D^*\ along with the character- 
istics of the two-body resonance, may be reflected in the distribution in space and energy 
of the detected particles. Typically, (e.g., Ref. [s^) there will be an enhancement or de- 
enhancement of the total cross section relative to the situation where three body scattering 
does not occur in the final state. In particular, the behavior of the X(3872) produced in iso- 
lation should be distinguishable from its behavior when in the presence of 5* = D^, D*^, D^, 
or Z)*°, a situation that may be accessible in the rich environment of the LHC. For example, 
if to leading order we can assume that the X and D^*^ are produced isotropically in their 
opening angle, we might attribute violations of isotropy at small opening angles to the effect 
of the large XD^*^ scattering in the final state. 
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